Giant Isotope Effect in the Incoherent Tunneling Specific Heat of the Molecular 

Nanomagnet Fe 8 
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Time-dependent specific heat experiments on the molecular nanomagnet Fes and the isotopic 
enriched analogue 57 Feg are presented. The inclusion of the 57 Fe nuclear spins leads to a huge 
enhancement of the specific heat below 1 K, ascribed to a strong increase in the spin-lattice relaxation 
rate V arising from incoherent, nuclear-spin-mediated magnetic quantum tunneling (MQT) in the 
ground-doublet. Since V is found comparable to the expected tunneling rate, the MQT process has 
to be inelastic. A model for the coupling of the tunneling spins to the lattice is presented. Under 
transverse field, a crossover from nuclear-spin-mediated to phonon-induced tunneling is observed. 

PACS numbers: 75.40.-s, 75.45.+J, 75.50.Xx 



Single-molecule magnets are fascinating nano-size su- 
perparamagnetic particles which at low temperatures 
may flip their magnetic moments by magnetic quantum 
tunneling (MQT) through the anisotropy barrier pj . Ob- 
servation of quantum tunneling in these molecules illus- 
trates the complexity of the interaction of such magnetic 
qubits with their "environment" (neighboring particles, 
nuclear spins, phonons). Indeed, the tunnel splitting A of 
the magnetic ground-state is many orders of magnitude 
smaller than the energy bias £ from, e.g., dipolar interac- 
tions between molecules, rendering MQT impossible at 
first sight. It is by now well established, both theoreti- 
cally fli and experimentally ggQ, that incoherent 
MQT is yet possible through the presence of rapidly fluc- 
tuating nuclear spins. The resulting dynamical hypcrfinc 
bias may, at any time, bring a fraction of the molecu- 
lar spins into resonance, thus opening an energy window 
E w 3> A for incoherent tunneling. 

It should be emphasized that in the Prokof'ev/Stamp 
(PS) model 0, relaxation of the magnetic moment of 
the molecular cluster (hereafter, cluster spin) is to the 
nuclear spin bath. A coupling to the lattice is not con- 
sidered, the argument being pj that only at long times 
relaxation by phonons will become more efficient than 
the nuclear-spin-mediated magnetic relaxation. How- 
ever, time-dependent specific heat experiments 0, Q 
neatly show that (for high enough MQT rate) both the 
electronic and nuclear spin systems can remain in thermal 
equilibrium with the lattice even deep into the quantum 
regime, where the only fluctuations possible are those 
arising from MQT events. This strongly suggests that, 
whereas the nuclear spins are needed to relax the (other- 
wise blocked) electron spins through MQT, at the same 
time the MQT mechanism apparently enables relaxation 
of both nuclear and electron spins to the lattice. Inter- 
estingly, in magnetic insulating compounds such as these, 



relaxation of nuclear spins to the lattice has to occur via 
the electron spin-phonon channel, direct nuclear relax- 
ation to the lattice being extremely slow at low temper- 
atures. In other words, by enabling the cluster spins to 
tunnel, the nuclei themselves can relax to the lattice! 

In this Letter we present definite proof for this unusual 
scenario. By comparing the low-T specific heat arising 
from incoherent tunneling in the ground-doublet for Fes 
and its 57 Fe-enriched counterpart 57 Fes, we show that the 
inclusion of the 57 Fe nuclear moments in the (otherwise 
identical) molecules leads to an enormous enhancement 
of the specific heat below 1 K, which can only result 
from a strong increase in the spin-lattice relaxation rate 
r. Below 1 K, this rate is found to be several orders of 
magnitude larger than predicted for conventional spin- 
lattice relaxation [t| . To explain our results an extension 
of the PS model is presented that includes inelastic tun- 
neling events, in which a spin- flip is accompanied by the 
creation or annihilation of low-energy phonons. 

Low-temperature (0.1 K < T < 7 K) specific heat C 
measurements were performed in a home-made calorime- 
ter using the thermal relaxation method. By varying 
the thermal resistance of the link between calorimeter 
and cold-sink, the characteristic timescale r e of the ex- 
periment can be varied. In this way measurements of the 
time-dependent C can be exploited to probe the long- 
time (0.1 — 10 2 s) magnetic relaxation Q. At higher 
temperatures (T > 2 K), C was measured using a com- 
mercial calorimeter. The enriched sample, denoted by 
57 Feg, was prepared as described in Ref. 5]. Except for 
the presence of the 57 Fe nuclear moments in 57 Fes (en- 
riched to 95% in 57 Fe), no difference in the magnetic 
structure between Fes and 57 Fes is to be expected. 

We first discuss the zero-field specific heat data. In 
Fig. 1, we compare the specific heat C/R of a non- 
oriented powder 57 Fes sample, as measured with two dif- 
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FIG. 1: (color online) Zero-field specific heats of non-oriented 
samples of standard Fes and of 57 Fes as a function of temper- 
ature. For 57 Fe8, data for r e « 35 s (•) and 1 s (o) are given, 
whereas for Feg, r e « 0.5 s (*). Drawn curves are explained 
in the text. 



ferent thermal links, with our previous data [6( on non- 
oriented Fes- The specific heat contains contributions 
from the phonon modes of the crystal lattice as well as 
magnetic contributions from electronic and nuclear spins. 
The lattice specific heat can be described by the sum of a 
Debye curve, which describes the contribution of acous- 
tic phonon modes, plus an Einstein oscillator term that 
probably arises from intramolecular vibrational modes 
(optical phonons). The overall fit yields for both sam- 
ples D ~ 19 K and 0£ ~ 38 K for the Debye and Ein- 
stein temperatures, respectively [13 , (plotted in Fig. 1 
as dotted curves). 

Above 1 K, the magnetic contribution C m to the spe- 
cific heat of both compounds is independent of r e . In 
this range, the equilibrium C m mainly arises from tran- 
sitions between the energy levels of the molecular spin 
S = 10, split by the uniaxial anisotropy 0,0. The as- 
sociated multilevel Schottky anomaly Co is shown as the 
dash/dotted curve in Fig. 1. Below 1 K, the remaining 
entropy arises almost completely from the ground-state 
doublet ±10, which is split by intercluster magnetic dipo- 
lar couplings and by hyperfine interactions between the 
electronic and nuclear spins. Dipolar interactions can 
induce, under equilibrium conditions, a long-range or- 
dered magnetic phase. In Fig. 1, we plot the equilibrium 
specific heat C eq of Fes, obtained from Monte-Carlo cal- 
culations 0| ! which predicts the occurrence of such a 
phase transition at Tc — 0.43 K. However, the equilibra- 
tion of the relative populations of the two lowest levels, 
either by thermal activation or by tunneling, is a very 
slow process 0] , thus C m measured at finite r e does not 
reach its equilibrium value. Indeed, as shown in Fig. 1 
(see also Q), the specific heat deviates from equilibrium 
below a blocking temperature Tb — 1.2 K, indicating 



that the electronic spin-lattice relaxation rate T becomes 
smaller than t~ . 

In this non-equilibrium regime we observe a spectacu- 
lar isotope effect. For the standard Fes, C m ~ Cq below 
T ~ 1 K, decreasing exponentially to nearly zero (be- 
low ~ 0.4 K the measured C basically equals the Debye 
term). By contrast, C m of 57 Fes is about 100 times larger 
and it increases with r e . We mention here that r e varies 
weakly with T and is determined for each data point. 
The T e values mentioned in Figs. 1 and 2 are those es- 
timated at T = 0.2 K for each particular T-sweep. The 
strong isotope effect is a direct evidence that T, i.e. the 
thermal contact between the lattice and electronic spin 
systems, is enhanced as a result of the introduction of 
57 Fe nuclear spins. Although the rate is not yet suffi- 
cient to ensure complete thermal equilibrium within the 
experimental time constants available, a sizable part of 
the entropy of the ground-doublet is now removed. 

It should be added that the 57 Fe nuclear spins 
1 = 1/2 will also contribute to the zero-field C/R. 
In thermal equilibrium, this contribution amounts to 
C hf { 57 Fe)/R = A 2 s 2 I(I + l)/3(fc s T) 2 , where A is the 
hyperfine coupling and s = 5/2 is the Fe 3+ electronic 
spin. Taking A/ks = 1-65 mK estimated by Stamp and 
Tupitsyn [lj|], we obtain the T _2 -term indicated by the 
lower dashed curve in Fig. 1. The higher dashed curve 
in Fig. 1 gives the equilibrium nuclear Chf calculated by 
adding the contributions of nuclear spins at the 120 pro- 
tons, the 18 14 N nuclei and the 8 79 > 81 Br nuclei present 
in the Fes molecule. Note that at T ~ 0.2 K, Chf is one 
to two orders of magnitude smaller than the measured 
C for 57 Fes. Further, see for standard Fes that C be- 
comes smaller than Chf at T w 0.1 K, suggesting that 
also nuclear spins are off equilibrium. 

We next discuss experiments performed under applied 
magnetic field. Figure 2 shows results for 57 Fes mea- 
sured at T ~ 0.22 K, compared to previous data for 
Fes obtained for nearly the same T and r e 6]. At 
such high-B and low-T, C m of a randomly oriented sam- 
ple is dominated by the contribution of those crystals 
whose anisotropy axes lie nearly perpendicular to the 
field. To substantiate this statement, we calculated that, 
at T = 0.2 K and B = 1.5 T, those crystals making 
an angle smaller than 87.5° with the field contribute less 
than one percent to the electronic equilibrium C. These 
experiments give thus information on the way the spin- 
lattice relaxation is modified by B± . 

In agreement with the zero-field-behavior, the C m 
curves of the two isotopic derivatives are very different 
for B < 1.5 T (Fig. 2). For higher fields, however, they 
are seen to merge within the experimental uncertainties. 
The observed dependence on the applied field can be well 
explained in terms of a tunable quantum tunneling rate. 
The perpendicular field introduces off-diagonal terms in 
the spin Hamiltonian, which increase the tunnel split- 
ting A by many orders of magnitude. When A becomes 
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FIG. 2: (color online) Field-dependence of the specific heat 
of 57 Fe 8 measured at T = 0.22 K, (for r e ^ 35 s), together 
with previous data |(J obtained for Fes at 0.24 K. Inset: Cal- 
culated T(B±) associated with 'conventional' direct processes 
(see text) for several T. 



FIG. 3: (color online) Spin- lattice relaxation rates of standard 
Fes and 57 Fes obtained by different experimental techniques, 
as labelled. Dashed curve is the Arrhenius fit to the high-T 
data. Arrows indicate low-T limits deduced from magnetiza- 
tion [j| and NMR data. 



larger than E w , the tunneling rate becomes no longer de- 
termined by the hyperfine interactions but rather by A. 
Then T should become nearly the same for both com- 
pounds, as observed. Besides phonon-assisted tunneling, 
the increase of A also leads to strong enhancement of 
the relaxation rate associated with 'conventional' direct 
processes of emission and absorption of phonons. This is 
shown in the inset of Fig. 2, where we plot T(B±) calcu- 
lated 14] using the anisotropy parameters from Ref . 15] . 
For B± > 2.2 T, T becomes of the same order of the 
experimental r~ . This explains why conventional the- 
ory for spin-lattice relaxation accounts quantitatively for 
the transition to equilibrium observed at large transverse 
fields |f], although it completely fails to account for the 
zero field data. 

To obtain quantitative information on the rate at 
which the spins approach thermal equilibrium at low T 
in zero field, we need to assume a particular expression 
for the time-dependence of C m . This is complicated by 
the fact that, at low T and especially below the order- 
ing temperature Tq, relaxation to equilibrium becomes 
a collective process in which each spin-flip modifies the 
dipolar biases acting on the other spins 0- This prob- 
lem was theoretically studied by Fernandez 0], who 
calculated numerically the time-dependent C m of a lat- 
tice of interacting Ising spins flipping by quantum tun- 
neling. Within this model, C m ~ C + (c T%/T)(T/v) 
where v = dT/dt is the temperature sweeping-rate and 
c is a constant that depends on the symmetry and lat- 
tice parameters. In our experiments, v ~ AT/r e , where 
AT ~ 0.05 T is the T-change of the calorimeter in each 
data point. It is therefore possible to determine T, up 
to a constant factor, from C m . Moreover, above Tc 
we can also fit C m assuming a simple exponential de- 



cay, C m (t) = C eq + (C - C eq ) exp(-ir) putting t = r e , 
where C eq is the calculated equilibrium C of the elec- 
tron spins (see Refs. 0, 0] for details). As shown in 
Fig. 3, the rates obtained by these two methods overlap, 
giving c = 0.05 ks- Deducing T in this manner is how- 
ever restricted to T regions for which the measured C m 
is sufficiently large compared to Cq and to the nuclear 
contributions f . Unfortunately, this is not the case for 
standard Fe 8 below 1 K. 

Figure 3 shows T(T) estimated as above together with 
data obtained at higher temperatures from Ref. 0]. 
Above 1 K, r follows an Arrhenius law that is approx- 
imately the same for both Fes derivatives and corre- 
sponds to a thermally activated relaxation over a barrier 
U w 22.5 K. By contrast for T < 1 K, the decrease of V of 
57 Fe8 abruptly slows down and thus T becomes orders of 
magnitude faster than predicted for activated behavior. 
Although we cannot extract T directly from the low-T C 
of standard Fes, we can still estimate an upper bound 
for it. Using the experimental r e (T) of standard Fe§ 
and assuming T(T) of Fe 8 to be directly proportional to 
r(T) of 57 Fe 8 , we calculate C m (T) with either the model 
of Ref. or exponential decay (see above). The so- 
obtained C m (T) reproduces well the experimental C m (T) 
of standard Fes if r(T) is taken a factor of 3 smaller than 
that of 57 Fe 8 . The same factor was obtained from time- 
dependent magnetization experiments ||, which pro- 
vided as well a similar low-T limit of T, also plotted in 
Fig. 3. These data have been successfully interpreted 
in terms of the PS model, i.e. nuclear spin-mediated 
tunneling events at a tunneling rate l/r t , followed by 
square-root relaxation through redistribution of dipolar 
fields throughout the sample in combination with addi- 
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tional spin-flips. In the PS model the tunneling rate is 
given by l/r t » A*/E w . With A t w 5 x 1CT 8 K and tun- 
neling window E w « 0.03 K, as found for 57 Fe 8 |{|, one 
obtains l/r t « 1 x 10~ 3 s _1 . The square root relaxation 
rate is basically given by [U 1/tq ~ (l/T t )(E w /E dip ) 2 . 
With E w « 0.03 K and E dip w 0.1 K, it follows that 
tq /r t « 10, in agreement with the experimental mag- 
netisation relaxation rate of 10~ 4 s _1 observed at low- 
est T El . As argued by Morello et al. 0] and Baek et 
al. |l9| , the quantum tunneling fluctuations of the clus- 
ter spins should set a T-independent lower limit to the 
longitudinal nuclear relaxation rate 1/T™ equal to the 
quantum tunneling rate. Indeed, the calculated value of 
l/r t is precisely that attained below about 0.5 K by the 
the experimental l/T? for 57 Fe 8 19] (Fig. 3). 

Our specific heat data thus neatly confirm the isotope 
effect predicted by the nuclear spin-mediated tunneling 
(PS) model and previously seen in Ref. 0, but in ad- 
dition provide strong evidence that in these incoherent 
quantum tunneling processes an efficient coupling to the 
phonons has to be involved, i.e. relaxation of the cluster 
spins is actually to the lattice at a rate corresponding 
to the square root relaxation, and not towards the nu- 
clear spin-bath, as assumed in the PS model. However, 
the problem is that, due to the strong average dipo- 
lar bias Edip ~ 0.1 K, combined with the very small 
l/r t « 1 x 10~ 3 s _1 , a direct coupling of the tunneling 
levels to the phonons via the usual spin-lattice interac- 
tion by phonon modulation of the crystal field leads to 
astronomically long relaxation times (inset of Fig. 2). We 
thus propose a two-step relaxation process as an alterna- 
tive. All data show that in zero field, in order to flip 
the spin, we have to rely on the dynamic hyperfinc inter- 
action. Through intercluster nuclear spin diffusion [lif 
the hyperfine bias can fluctuate over an appreciable part 
of the hyperfine split manifold. Since the width of this 
energy window E w is of the same order (0.01 to 0.1 K) 
as the dipolar bias, near-resonant conditions for the tun- 
neling levels can be met by the combined action of both 
biases. Importantly, (i) when the spin flips, both the 
hyperfine field and the dipolar bias acting on it change 
sign, implying that in this process energy can be inter- 
changed between nuclear spins and electron dipolar inter- 
action reservoir, and (ii) the reversal of this dipolar bias 
is instantaneous (picosec) compared to all relevant time 
scales and produces a temporary disturbance in the local 
dipolar field distribution. But since the cluster spins are 
on the nodes of a crystal lattice, they are coupled not 
only by dipolar but also by weak intercluster elastic (van 
der Waals) forces. Following general arguments on en- 
ergy and angular momentum conservation |20l | , this may 
result in a temporary local lattice instability, similar as 
in a Franck-Condon type electronic transition associated 
with light absorption by lattice defect or impurity states 
in solids [2l| . Thus an unbalance of the interchanged hy- 
perfine and dipolar energy quanta can be taken up by the 



lattice as potential energy in the form of a local phonon 
mode, followed by dissipation into thermal phonons and 
simultaneous outward evolution of the dipolar field dis- 
tribution from the 'defect'. The thermalisation of local 
vibrational modes due to anharmonic processes has been 
studied by several authors 0, ^ and the associated 
times are estimated to be of order 10 3 to 10 4 s at the low 
temperatures (« 0.1 K) considered here [13], faster in- 
deed than the observed spin-lattice relaxation times (de- 
termined by the tunneling rate). 
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